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ABSTRACT 

Most Galaxy-sized systems (Mh os t — 10 12 Af©) in the ACDM cosmology are expected to have in- 
teracted with at least one satellite with a total mass Af sat ~ 10 11 M Q ~ 3Mdi s k in the past 8 Gyr. 
Analytic and numerical investigations suggest that this is the most precarious type of accretion for the 
survival of thin galactic disks because more massive accretion events are relatively rare and less mas- 
sive ones preserve thin disk components. We use high-resolution, dissipationless iV-body simulations 
to study the response of an initially-thin, fully-formed Milky- Way type stellar disk to these cosmolog- 
ically common satellite accretion events and show that the thin disk does not survive. Regardless of 
orbital configuration, the impacts transform the disks into structures that are roughly three times as 
thick and more than twice as kinematically hot as the observed dominant thin disk component of the 
Milky Way. We conclude that if the Galactic thin disk is a representative case, then the presence of 
a stabilizing gas component is the only recourse for explaining the preponderance of disk galaxies in 
a ACDM universe; otherwise, the disk of the Milky Way must be uncommonly cold and thin for its 
luminosity, perhaps as a consequence of an unusually quiescent accretion history. 
Subject headings: Cosmology: theory — galaxies: formation — galaxies: evolution 



1. INTRODUCTION 

A solid majority of observed galaxies have disk- 
dominant morphology; despite wide variance in methods 
of sampling and classification, roughly 70% of Galaxy- 
sized dark matter halos in the universe host late- type sys- 
tems (IWeinmann et al.l 120061: Ivan den Bosch et al 1120071 : 
lllbert et all 120061: lChoieTal.1 120071 : iPark et all 120071 
Moreover, iKautsch et al.l ( 2006( 1 find that about one- 
third of all local disk galaxies have no observable 
pressure-supported component (whether a "classical" 
bulge formed by the central starburst associated with 
a merger event, or a "pseudobulge" having arisen from 
the secular transport of angular momentum towards the 
galactic center), and another one-third host systems with 
only pseudobulges, a conclusion supported by spheroid- 
disk decomposition of larg e galaxy samples (|Allen et al.1 
2006; Ba razza et al.l 120081 ) . The vast majority of disk 
stars in the Milky Way reside in the thin disk compo- 
nent, with an exp onential scale height of za — 300 ±60 pc 
(|Juric et al.ll2008l and references therei n) and a total ve- 
locity dispersion of a to t — 35 km s _1 (|Nordstrom et al.1 
2004). Whether the scale height of the Galactic disk is 
typical for galaxies of its size is a topic of vital inter- 
est. Unfortunately, firm measurements for a statistical 
sample of galaxies have been limited by dust obscura- 
tion, which present a problem even in K-band imaging 
(|Kreeel et alj|2005t lYoachim fc Dalcantonll2006[) . 

Aside from the considerable challenges associated 
with forming disk galaxies in ACDM cosmologies (e.g., 
iMaver et al.ll2008l ). hierarchical models must also self- 
consistently maintain thin, rotationally-supported sys- 
tems against the constant barrage of merging subhalos. 
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Though the former end e avor h as enjoyed some recent 
advan ce s (lAbadi et al.1 12003 : ISommer-Larsen et al.l 
20031: iBrook et all l200l iRobertson et all 12004 
Governato et al. I 120071) . the survival of disk galax- 



ies during the often-violent mass acc retion history of 
their dark host remains a concern (iToth fc Ostriker 
1991 IQuinn et all Il99l iWalker et all 119961 : lWv~ 
200 lh and has been the target of numerous studies 
aimed at quantifying the r esilience of galactic disks 
to satellite accretion e vents (IQuinn fc Goodman) [1 986: 
Toth fc Ostriker! [T992: Quinn et alj 119931: IWalker et al 
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Both numerical simulatio ns (IStewart et al.l 120081) and 
purel y analytic calculations (jPurcell et al.ll2007l : IZentnerl 
I2007h indicate that mass delivery into dark matter halos 
of mass Mhost is dominated by th e accretion of o bjects 
with mass - (0.05 - 0.15)M host . IStewart et ail (|2008l ) 
find that ~ 70% of 10 12 M© Galaxy-sized halos have ac- 
creted a system of mass M sat ~ 10 11 M© ~ 3Mdi s k into 
their virial radii in the last 10 Gyr, with associated disk 
impacts within the last 8 Gyr. IStewart et al.l (2008) also 
find that less massive accretions are virtually ubiquitous, 
and that the merger fraction falls off quickly for satellites 
larger than M sat ~ 2 x lO 11 ^/©. Overall, these results 
suggest that ~ 1 : 10 satellite accretion events represent 
the primary concern for disk surv ival in ACDM. Along 
these lines, recent analytic work bv lHopkins et al.1 (2008) 
suggests that orbital energy deposition via merger is less 
destructive to a disk than was often previously surmised 
(IToth fc Ostriker! fl99l [Quinn et al.lll993l: IWalker et all 
1996), claiming that the Galaxy could have undergone 
^5 — 10 independent mergers of this kind since z ~ 2 

while main taining a thin disk . 

Recently, Ka zantzidis et al.l (|2008f) utilized dissipation- 
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Fig. 1. — Edge-on surface brightness maps, assuming M+/L = 3, for primary galaxies 1 (upper panels) and 2 (lower panels). Initial 
models (t = Gyr) are shown in the left panel, while the results (t = 5 Gyr) for satellite-infall orbital inclinations of 9 = 30° and 90° 
appear in the center and right panels, respectively. 



less iV-body simulations to investigate the response of 
thin galactic disks subject to a ACDM-motivated satel- 
lite accretion history. These authors showed that the 
thin disk component survives, though it is strongly per- 
turbed by the viole nt gravitational encounte rs with sub- 
structure. However. iKazantzidis et all (|2008l ) focused on 
infalling systems with masses in the range 0.2Mdi s k ^ 
-Wsat ^$ A^disk, ignoring the most massive accretion events 
expected over a galaxy's lifetime. In this Letter, we ex- 
pand upon this initiative by investigating the morpho- 
logical and dynamical evolution of initially-thin Galaxy- 
type disks during cosmologically common ~ 1 : 10 ac- 
cretion events involving two-component (stars and dark 
matter) satellites of mass M sat ~ 10 n M w ~ 3Mdj s k- 

Wor king in a similar mass regime, IVillalobos fcH elmi 
(2008) simulated the formation of thick disks via the in- 
fall of satellite galaxies with virial masses ~ 10—20% that 
of the host, using both az = Galactic primary system 
and a scaled version at z — 1 in order to show that real- 
istic thick disks result from these impacts. Though our 
preparation is similar, our goals and techniques are differ- 
ent. We aim to determine whether any thin, dynamically 
cold component can survive such an event, and conserva- 
tively use a primary disk that is as massive as the Milky 
Way disk today. 

Past investigations into the stability of galactic 
disks against the infall of satellites have often suf- 
fered from the necessities of numerical limitations 
or from analytic axioms later deemed incompati- 
ble with standard cosmological models; for exam- 
ple, the modeling of one or more structural com- 
ponents as r igid p o tentials (iQuiim fcOoodmanl HoM 
Quiim etaLl 119931: ISellwood et a ll 119981: l Ardi et all 
20031 lHavashi fc Chibal I2006D . the initialization of 
a disk much thicke r than the old, thin stellar 
disk of the Galaxy dQuinn et al.l [l993u IWalker et al.1 
199llHuang fc Carlberdll997tlVelazauez fc Whitd ll999l: 
Font et alJ 120011 : IVillalobos fc Helmil l2008f) . "the in- 
fall of satellite s with only a concentrated bary- 
onic component Jo umn "et^Ll[l99l IWalker et al. | [ l996: 
iHuang fc Carlberd 119971 : IVelazquez fc White! Il999f) . 
and the imposition of subhalo infalls with or- 
bital parameters inconsistent with ACDM cosmolog- 



ical models dQuinn et al.1 119931 : IWalker et afl 119961 : 
IHuang fc Carlberd Il997fh Analytic arguments, mean- 
while, have historically been forced to assume simplifica- 
tions su ch as the local deposit ion of a satellite's orbital 
energy dToth fc Ostrikerlll992D , or the absence of globa l 
heating modes (jBenson et al.ll2004t iHopkins et al.ll2008h 
which are analy t ically shown to dominate disk heating by 
ISellwood et al.l (|1998l ). although the latter authors em- 
ploy a rigid satellite model and perfectly radial polar or- 
bits for their simulated experimental tests. Fortunately, 
advances both in computational power and in our un- 
derstanding of ACDM expectations allow us to address 
these concerns directly. 

Our contribution improves upon earlier studies in sev- 
eral important respects. First and foremost, we examine 
the response of galactic disks to accretion events that rep- 
resent the primary concern for disk survival in ACDM 
cosmologies. Secondly, we employ galaxy and satellite 
models that are constructed in equilibrium from fully 
self-consistent distribution functions and which have the 
resolution in force and mass to study the heating of a disk 
that is as thin as the old thin stellar disk of the Milky 
Way (zd ~ 300 pc); in synergy with the high mass and 
force resolution we adopt, this quality allows us to con- 
struct equilibrium Af-body models of disk galaxies that 
are as thin as the old, thin stellar disk of the Galaxy. 
Lastly, the masses, density structure, stellar content, and 
orbital configurations of our infalling satellites are di- 
rectly motivated by the prevailing ACDM paradigm of 
structure formation. 

2. METHODS 

All simulations are performed using t he multi - step ping, 
parallel, tree A^-body code PKDGRAV [Stadejl (|2001f L in 
which we set the gravitational softening length to e = 
100 pc and 50 pc for dark matter and stellar particles, 
respectively. 

2.1. Primary and Satellite Galaxy Models 

We construct A-body realizations of primary disk 
galaxi es and satellites using the method of Widrow et al. 
(2008). This technique produces self-consistent, multi- 
component galaxy models that are ideal for studying 
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complex dynamical processes associated with the intrin- 
sic fragility of galactic disks such as gravitational inter- 
actions with infalling subhalos. We explore two initial 
models for the primary galaxy in our satellite-disk en- 
counter simulations: Galaxy 1 (hereafter Gl), a Milky- 
Way-analog system drawn from the set of self-consistent 
equilibrium models that bes t fit Galactic observa tional 
parameters as produced by IWidrow et alj (|2008f ): and 
Galaxy 2 (hereafter G2), an identical system save for 
the absence of a central bulge, i.e., the two models have 
stellar disks and dark halos with equivalent initial prop- 
erties. In each case the dark matter halo of the pri- 
mary gala xy was popula t ed by 4 x 10 6 particles fol- 
lowing the iNavarro et alj (|1996l hereafter NFW) den- 
sity profile with scale radius r s — 14.4 kpc, and the 
bulge in Gl (comprised of 5 x 10 5 particles) contained 
a stellar mass Mbuige = 9.5 x 10 9 M Q following a Sersic 
profile with effective radius R e = 0.58 kpc and index 
n = 1.118. The stellar disks, comprised of 10 6 par- 
ticles each, contained a mass A/disk = 3.6 x 10 10 M Q 
following an exponential distribution in cylindrical ra- 
dius with scale length Rd — 2.84 kpc, while the vertical 
distribution of stars was described by a sech 2 function 
with Zd = 0.43 kpc being the vertical scale height. We 
note that the choice of numerical and physical parame- 
ters minimize secular evolution (e.g., strong bar forma- 
tion, artificial heating through interactions with massive 
halo particles) on the timescales of relevance to our in- 
vestigation, which could interfere with the interpretation 
of our results. In the left panel of Figure [IJ we show 
the edge-on surface brightness map for both primary 
galaxy models, having assumed a stellar mass-to-light 
ratio M+/L = 3. The satellite galaxy in each case was 
initialized with 9 x 10 5 dark particles representing a mass 
Msat = 1.0 x 10 11 M Q within the virial radius of a halo 
which is well-fit by an NFW profile with a concentration 
c v i r ~ 14 at z — 0.5. We populate this satellite with 
a stellar mass M± = 2.2 x 1O 9 M , roughly correspond- 
ing to the upper- la limit derived bv lConrov fe Wechsleri 
(2008) for M+/M Bat (z ~ 0.5) at our subhalo's virial mass, 
and we distribute these 10 5 stellar particles in a cen- 
tral spheroid with Sersic index n ~ 0.5 according to the 
distribution of shape para meters versus dwarf e lliptical 
galaxy magnitudes found bv lvan Zee et all (|2004h in their 
survey of Virgo cluster members. 

2.2. Satellite Galaxy Orbits 

Our initial subhalo velocity vectors are motivated by 
cosmological investigations of substructure accretions, 
where the distributions of radial and tangential veloc- 
ity components (v r and v t ) peak respectively at 90% 
and 60% of the y i rial velocity of the satell ite's host 
halo ()Benson|[2005t iKhochfar fe Burkertj|2006f) . In our 
case this corresponds to an initial subhalo velocity with 
v r = 116 km/s and v t — 77 km/s. We initiate the infall 
of each simulation's subhalo at a relatively large radius 
of approximately 120 kpc to ensure that the disk does 
not suffer substantial perturbations due to the sudden 
presence of the satellite's potential well. We simulate an 
array of orbital inclinations (9 = 0°, 30°, 60°, and 90°, 
defining 9 as the angle between the angular momentum 
axes of the disk and the orbit) in order to assess the con- 
sequence of this parameter on the evolution of the galac- 
tic disk. In the polar infall (9 = 90°), we eliminate the 
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Fig. 2. — Minor-axis surface brightness profiles for initial 
final models at two Galactocentric radii: R = Rq = 8 kpc 
panel) and R = 2Rq = 16 kpc (right panel). 
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Fig. 3. — The thin- and thick-disk scale heights in the final state 
(t = 5 Gyr) for e ach of our s imu lated galaxies, compared to the 
values derived bv lJuric et alj 1120081 ) for the Milky Way. The two 
panels show the result of a two-component sech 2 fit, with the upper 
(lower) panel describing the thin (thick) disk's scale height. 

tangential velocity component of the subhalo, sending 
the satellite on a direct-impact trajectory into the center 
of the primary galaxy; this case is somewhat unrealistic, 
but provides an interesting experimental benchmark. All 
but one of the non-polar subhalo orbits are initialized as 
prograde with respect to the primary galaxy's rotation; 
we also simulate a retrograde orbit for Gl with 9 = 60° 
in order to investigate wh ether the heating effects a re re- 
duced (as conjectured bv lVelazquez fc Whitdll999f ). All 
simulations were evolved for a total of 5 Gyr, after which 
the subhalo has fully coalesced into the center of the 
host halo and the stellar disk has relaxed into stabil- 
ity; although there are certainly remnant features in the 
outer disk and halo that will continue to phase-mix and 
virialize on a much longer timescale, our investigations 
indicate that the disk-evolution process has reached a 
quasi-steady state by this point in the encounter's evo- 
lution. 

3. RESULTS AND IMPLICATIONS 

Edge-on surface brightness profiles for remnants of the 
9 = 30° and 90° impacts are shown in the middle and 
left panels of Figure 1 , where the upper and lower render- 
ings correspond to primary cases Gl and G2, respectively 
(with and without initial bulge). It is clear from these 
images that the resultant disks are considerably thicker 
than the initial case. We note that while the stars in 
the accreted satellite end up in the final disk remnant 
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Table 1 

Final (t = 5 Gyr) Galaxy Properties at R = Rq = 8 kpc 



Orbital 


2 thin 


z thick 


< 1*1 > 


-^median 


Rd 


Cz 


Ctot 


Inclination 


(Gl, G2; 


(Gl, G2; 


(Gl, G2; 


(Gl, G2; 


(Gl, G2; 


(Gl, G2; 


(Gl, G2; 


of Subhalo 


in kpc) 


in kpc) 


in kpc 


in kpc) 


in kpc) 


in km/s) 


in km/s) 


initial (t = Gyr) 


0.43,0.43 


N/A 


0.3,0.3 


0.3,0.3 


3.0,3.0 


19.1, 18.7 


50.9,52.1 


9 = 0° (prograde) 


0.80,0.97 


3.70,3.65 


0.9, 1.0 


0.5,0.6 


2.3,4.5 


25.1,28.0 


115.2, 107.1 


= 30° (prograde) 


1.68, 1.75 


5.10,5.30 


1.7,1.8 


1.0,1.1 


3.5,2.9 


37.9,40.6 


95.5,102.6 


9 = 60° (prograde) 


1.33, 1.30 


4.38,4.35 


1.8,2.0 


0.9, 1.1 


2.2,2.6 


33.5,35.1 


82.4, 86.1 


e = 60°-retro (Gl only) 


1.18 


6.50 


2.1 


0.8 


2.6 


31.4 


83.3 


8 = 90° (polar) 


1.05, 1.08 


9.35,9.45 


2.0,1.9 


0.6,0.7 


4.2,3.0 


26.2,29.4 


70.0,75.6 


Milky Way (observed)" 


0.34 


1.01 


0.298 


0.208 


2.6,3.6 


~ 10 - 20 


~ 30 - 40 



(thin, thick) 

"For the Galaxy's em pirical constraints, we quote the disk scale heights and lengths derived by Uuric et al.l (120081) and the velocity 
dispersions obtained by [Nordstrom ct al. (2004) for solar- neighbor hood stars of median age (t ~ 2 — 3 Gyr). 



(c.f. Villalobos & Hclmi 2008), primary disk stars domi- 
nate these images, even high above the plane. 

Figure [5] shows the minor-axis surface brightness pro- 
files for the Gl simulations using = 3. The left 
panel shows a vertical slice at a projected radius of 
Rq = 8 kpc and the right panel shows a similar slice 
at radius 2 Rq . Black solid lines show the initial disk 
and different color/line types represent the remnants as 
indicated. Clearly, the resultant disks are dramatically 
thicker than the initial galaxy model in each case. In 
order to conservatively compare our disks to the Milky 
Way, we allow for thick and thin components by fitting a 
double-sech 2 profile at Rq. The fitted scale heights are 
shown in Table[T]and comp ared directly in Fig ure [3] to the 
Galactic values obtained bv lJuric eFaDdQQi) 3 . Though 
we initially employ a disk that is thicker (z t hm = 0.43 
kpc) and therefore conservatively more robust to accre- 
tion events (Kazantzidis et al. 2009, in prepar ation) than 
the G alactic value of Zthin = 0.34 kpc from lJuric et"aT1 
(2008), the final systems all have thin-disk components 
with Zthin larger by a factor of ~ 3 — 5 than the Milky 
Way. Moreover, the low-surface-brightness thick com- 
ponent in our remnant disks is also considerably thicker 
than the Galactic thick disk, with scale heights so large 
(zthick — 4 — 10 kpc) that this material would likely be 
considered a stellar halo component. 

A second relevant measure of disk survival is the stel- 
lar velocity dispersion; we therefore compare the ve- 
locity ellipsoid of our final disks to that observed in 
the solar neighborhood by the Gen eva-Copenhagen Sur- 
vey ( Nordstrom et al.ll2004l see also lSeabroke fc Gilmord 
I200J). In Figure gl we show these values for veloc- 
ity dispersion in each coordinate, where the indicated 
range spans the stellar population age, and the point 
is placed at the median a ge t ~ 2 — 3 Gyr according 
to INordstrom et al.l (|2004f ) for this local sample. Shown 
also are the corresponding dispersion components for our 
initial and final stellar disks measured within an 0.3-kpc 



box centered on the disk plane at Rq = 8kpc. As sum- 
marized in Table 2, each of our simulated merger rem- 
nants are substantially enhanced in all three components 
of velocity dispersion {ur^^ corresponding to au,v,w)- 
The total dispersion cr tot = {a\ + er^ + tx 2 ) 1 / 2 increases 
by a factor of ~ 1.5 — 2 compared to that of the initial 
disk. 

4. CONCLUSIONS AND DISCUSSION 

Using fully self-consistent TV-body simulations of 
satellite-disk interactions we have quantitatively demon- 
strated for the first time that cosmologically common ac- 




MW obs initio! 0=0° 0=30° = 60° 0= 



90° 



3 The scale heights derived by Juric ct al. (2008) belonged to 
exponential profiles; we have therefore multiplied these values by 
a factor of 1.12 to obtain scale heights belonging to sech 2 profiles 
that fall by one mag/arcsec 2 at the same height as the exponential 
fits. This multiplicative factor is more appropriate for thin-disk 
comparisons near the peak of the profile than the widely- used factor 
of 2 that matches exponential and sech 2 profiles at large heights 
above the disk plane. 



Fig. 4. — The radial and vertical components of velocity dis- 
persion Ufi and a z (top and middle panels), as well as the to- 
tal stellar dynamical temperature otot, at the solar neighborhood 
(Rq = 8 kpc) of our simulated disks, compared to the local val- 
ue s obtained by the Gene va-Copenhagen survey results described 
in INordstrom et~a l. (2004). In each coordinate, the observational 
spread is marked by a dotted line and the dispersion of the sample's 
median-age stars (t ~ 2 — 3 Gyr) is denoted by a diamond. 
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cretion events of mass ratio ~ 1 : 10 do not preserve thin, 
dynamically cold stellar disks like the old, thin stellar 
disk of the Milky Way. This has potentially serious ram- 
ifications for models of galaxy formation and evolution. 
It is possible that our benchmark case of the Milky Way 
is not representative, and that the Galaxy sits within a 
rare halo that has not experienced in the last ~ 8 Gyr a 
disk impact associated with a significant accretion event, 
as posited in the obse rvationally-motivated suggestion of 
lHammer et al.l (|2007T ). in which the Galaxy is shown to 
have remarkably low angular momentum and stellar mass 
compared to local spiral galaxies in host halos of similar 
mass. Future investigations may help quantify the range 
of thin-disk scale heights in the local universe. 

Otherwise, the addition of gas physics may play a role 
in explaining the apparent discrepancy. Gas can cool and 
reform a thin d isk, and its presence ma y stabilize the stel- 
lar disk (e.g., iRobertson et alj [20061 ). The regrowth of 
the massive thin disk after a satellite accretion may cause 
heated stars to contract and lose kinetic energy. Accu- 
rate treatment of the various aspects of hydrodynamics 
will therefore play a crucial role in the capacity of sim- 
ulated galaxy evolution to reproduce thin disks such as 
those that dominate observed galax y catal ogs. 

In a recent paper iHopkins et alj (|2008f ) have argued 
that disk heating is less effective than previously thought 
and that the expected merger histories of ACDM halos 
are compatible with the high thin-disk fraction seen in 
the Universe. It is important, therefore, to investigate 
this point of disagreement. Their result, a reshap ing of 
the arguments presented in iToth &; Ostrikerl (|1992f ) (up- 
dated to reflect the more realistically radialized orbits of 
a ACDM cosmology), relied primarily on an analytic for- 
mula, normalized to simulations with much lower mass 
and force resolution than those explored here, to map the 
ratio (M sa t/Mhost) to a disk heating parameter AH/R, 
where H is the median scale height of the resultant disk 
and R is the radius where the height is measured (and 
must be within a factor of two of the disk half-mass ra- 
dius R e ). For the ~ 1 : 10 mass-ratio a ccretion events we 
explore here, the iHopkins et alJ (|2008l ) formula predicts 
a disk thickening of AH/R ~ 0.015. Our simulations 



typically exhibit significantly more heating; at R = R e 
we measure AH/R ~ (0.03 — 0.09) and, because of the 
impact-induced flaring, we measure even larger values 
AH/R ~ (0.05 - 0.11) at R = 2R e . 

It is perhaps not surprising that our results disagree 
with first-order analytic expectations. In addition to 
direct heating, the resultant disk structure is affected 
by global modes such as bending and de nsity waves ex- 
cited in the disk as the interaction occurs (Scllwoo d et alJ 
1998), and not included in the simple analytic scalings is 
a dependence on the orbital inclination of the encounter 
that is likely associated with resonant coupling. Finally, 
though IHopkins et alJ (|2008fl normalized their results to 
numerical simulations, those initial disks were signifi- 
cantly thicker than the Galactic-type disk we have simu- 
lated, and were therefore more robust to tidal perturba- 
tions. Direct numerical experiments involving satellite- 
disk encounters indicate that mass-ratio, orbital incli- 
nation, initial disk scale height, and relative dark matter 
fraction are all crucial in determining the degree to which 
galactic disks are perturbed by infalling subhalos (Fig- 
ures [3j |H and Kazantzidis et al. 2009, in preparation). 
More detailed analysis is forthcoming of the morpho- 
logical and dynamical effects experienced by our disks; 
among other concerns, we defer for future work the is- 
sues of stellar-halo/thick-disk distinguishability and the 
reinforcement of central bulges by accreted stars. 
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